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Abstract 

This work is devoted to the study of the initial boundary value problem for 
a general isothermal model of capillary fluids derived by J.E Dunn and J.Serrin 
(1985) (see [H]), which can be used as a phase transition model. We will prove the 
existence of strong solutions in finite time with discontinuous initial density, more 

precisely lnp is in B 2 2 00 (R JV ). Our analysis improves the results of 20i and 
[35] by working in space of infinite energy. In passing our result allow to consider 
initial data with discontinuous interfaces, whereas in all the literature the results of 
existence of strong solutions consider always initial density that are continuous. More 
precisely we investigate the existence of strong solution for Korteweg's system when 
we authorize jump in the pressure across some hypersurface. We obtain also a result 
of ill-posedness for Korteweg system and we derive a new blow-up criterion which is 
the main result of this paper. More precisely we show that if we control the vacuum 
(i.e - € L^(B^ +e 1 (R JV )) with e > ) then we can extend the strong solutions 
in finite time. It extends substantially previous results obtained for compressible 
equations. 

1 Introduction 

We are concerned with compressible fluids endowed with internal capillarity. The model 
we consider originates from the XlXth century work by Van der Waals and Korteweg 
[471, 137] and was actually derived in its modern form in the 1980s using the second gradient 
theory, see for instance [21, 361136]. The first investigations begin with the Young-Laplace 
theory which claims that the phases are separated by a hypersurface and that the jump 
in the pressure across the hypersurface is proportional to the curvature of the hypersur- 
face. The main difficulty consists in describing the location and the movement of the 
interfaces. 

Another major problem is to understand whether the interface behaves as a discontinu- 
ity in the state space (sharp interface) or whether the phase boundary corresponds to 
a more regular transition (diffuse interface, DI). The diffuse interface models have the 
advantage to consider only one set of equations in a single spatial domain (the density 
takes into account the different phases) which considerably simplifies the mathematical 
and numerical study (indeed in the case of sharp interfaces, we have to treat a problem 
with free boundary). 
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Another approach corresponds to determine equilibrium solutions which classically con- 
sists in the minimization of the free energy functional. Unfortunately this minimization 
problem has an infinity of solutions, and many of them are physically wrong (some de- 
tails are given later). In order to overcome this difficulty, Van der Waals in the XlX-th 
century was the first to add a term of capillarity to select the physically correct solu- 
tions, modulo the introduction of a diffuse interface. This theory is widely accepted as a 
thermodynamically consistent model for equilibria. 

Korteweg-type models are based on an extended version of nonequilibrium thermodynam- 
ics, which assumes that the energy of the fluid not only depends on standard variables 
but also on the gradient of the density. Alternatively, another way to penalize the high 
density variations consists in applying a zero order but non-local operator to the density 
gradient ( see [12], [33], [33] )■ For more results on non local Korteweg system, we refer 

also to pil EH EE]- 

Let us now consider a fluid of density p > 0, velocity field u G M. N , we are now interested 
in the following compressible capillary fluid model, which can be derived from a Cahn- 
Hilliard like free energy (see the pioneering work by J.- E. Dunn and J. Serrin in [21 j and 
also in [2j I1CH [2^]). The conservation of mass and of momentum write: 

( d 

—p + div(pu) = 0, 

f (1-1) 
— (pu) + div(pu ®u)- div(2p(p)D(u)) - V(A(p))divu) + VP(p) = divK, 

where the Korteweg tensor read as following: 

divK = V(pk(p)A p + ^(k( p ) + p K '(p))\Vp\ 2 ) - div ( K (p)Vp®Vp). (1.2) 

k is the coefficient of capillarity and is a regular function. The term div.ff allows to 
describe the variation of density at the interfaces between two phases, generally a mixture 
liquid-vapor. P is a general increasing pressure term. D(u) = (Vu+ Vit) being the stress 
tensor, p and A are the two Lame viscosity coefficients depending on the density p) and 
satisfying: 

p > and 2p + NX > 0. 

Here we want to investigate the existence of strong solution for the system (|l,ip when we 
authorize jump in the pressure across some hypersurface. To do this, we need to prove 
the existence of strong solutions in critical space for the scaling of the equations with 
initial densities which are not continuous. It will be one of main interest of this paper 
with a new blow-up criterion. 

Before entering in the heart of the subject, we now want to recall the classical energy 
inequalities. Let p > be a constant reference density (in the sequel we will assume that 
p = 1), and let II be defined by: 

so that P(s) = sll (s) — II(s) , II (p) = 0. Multiplying the equation of momentum 
conservation in the system (jl.ip by u and integrating by parts over R , we obtain the 
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following estimate: 



/ d p \uf + (n(p)-n(p)) + -L\Vp\ 2 )(t)dx+ 1 - f [ p( P )\D(u)\ 2 dxdt 

Jr" 2 2k P 2 Jo Jr n 3 ^ 

< / + (n(po) - n(p)) + -^|Vpol 2 )^, 

Jrn 2p 2kpq 
with: It follows that assuming that the initial total energy is finite: 

£o= [ (^ + (n(p )-n(p)) + ^|vpol 2 )dx<+oo, 

Jrn 2p 2 
then we have the a priori following bounds when P(p) = ap 1 : 

(p- 1) G L°°(Ll), and p\u\ 2 € L 1 (0, oo, L 1 (E 7V )), 

y^)VpG ^(0,00, L 2 ^))^, and G L 2 (0, oo, . 

We refer to section [2] for the definition of the Orlicz spaces. We now want to recall the 
notion of scaling for Korteweg's system. Such an approach is now classical for incom- 
pressible Navier-Stokes equation and yields local well-posedness (or global well-posedness 
for small data) in spaces with minimal regularity. Let us explain precisely the scaling of 
Korteweg's system. We can easily check that, if (p,u) solves (jl.ip . then (p\,U\) solves 
also this system: 

Px(t, x) = p(X 2 t, Xx) , u\(t, x) = Xu(X 2 t, Xx) 
provided the pressure laws P have been changed into X 2 P. 

Definition 1.1 We say that a functional space is critical with respect to the scaling of 
the equation if the associated norm is invariant under the transformation: 

(p,u) — > (px,u\) 

(up to a constant independent of X). 

This suggests us to choose initial data (po,uo) in spaces whose norm is invariant for all 
A > by (po,u ) — > Oo(A-), A«o(A-)). 



Strong solutions 

A natural candidate is the Besov space x {B^^ l ) N , but since B^^ is not included 

in L°°, we cannot expect to get L°° control on the density when po G B^ 2 . In particular 
it implies some difficulties in order ton control the non linear terms as the pressure. An 
other candidate would be (B^^ fl L°°) x (B^^ 1 ) N , however controlling the L°° norms 
of density is an hard task, even if the Korteweg system allows regularizing effects on the 
density. Indeed it appears also quadratic terms in the gradient of the density that it is 
necessary to control in L 2 (L°°) norm to be able to obtain estimate on p in L°°. 
This is the reason why in the literature the authors consider initial density which are 
in Banach spaces imbedded in Let us briefly mention that the existence of strong 
solutions for iV > 2 is known since the works by H. Hattori and D. Li [34, 35j. R. Danchin 
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and B. Desjardins in [20j improve this result by working in critical spaces for the scaling 

AT JV i 

of the equations, more precisely the initial data (po,poUo) belong to B 21 x B 21 (the 

JV 

fact that B 21 is embedded in L°° plays a crucial role to control the vacuum but also 
for some reason slinked with the multiplier space theory). In [38j, M. Kotschote showed 
the existence of strong solution for the isothermal model in bounded domain by using 
Dorea-Venni Theory and calculus. In [34J, we generalize the results of [20 in the case 
of non isothermal Korteweg system with physical coefficients depending on the density 
and the temperature. More precisely we get strong solutions with initial data belonging 

jv JV_ M. — 2 

to the critical spaces B 2l x B 21 x B 21 when the physical coefficients depend only 
on the density. 

Weak solutions 

We wish to prove existence and uniqueness results for (NHV) in functions spaces very 
close to energy spaces. In the non isothermal non capillary case and P{p) = ap" 1 , with 
a > and 7 > 1, P-L. Lions in [35] proved the global existence of variational solutions 
(p, u, 9) to (NHV) with k = for 7 > f if N > 4, 7 > ^ ^ N = 2,3 and initial data 
(po,mo) such that: 

n(po)-n(p), ^GL 1 ^), and e L 1 ^). 
Po 

These solutions are weak solutions in the classical sense for the equation of mass conser- 
vation and for the equation of the momentum. Notice that the main difficulty for proving 
Lions' theorem consists in exhibiting strong compactness properties of the density p in 
L v loc spaces required to pass to the limit in the pressure term P(p) = ap" ( . 
Let us mention that Feireisl in [23] generalized the result to 7 > y in establishing that we 
can obtain renormalized solution without imposing that p G Lf oc , for this he introduces 
the concept of oscillation defect measure evaluating the lost of compactness. 
We can finally cite the result from Bresch-Desjardins in [7], [8] where they show the exis- 
tence of global weak solution for (NHV) with k = and with a cold pressure. In fact 
they are working with specific viscosity coefficients which verify an algebraic relation. It 
allows them to get good estimate on the density by using new energy inequality and to 
treat by compactness all the delicate terms as the pressure. In [3D] . Mellet and Vasseur 
improve the results of Bresch,Desjardins by dealing with the case of a general pressure 
P(p) = ap" 1 with 7 > 1. 

In the case k > 0, we remark then that the density belongs to L°°(0, 00, ij^HX^)). Hence, 
in contrast to the non capillary case one can easily pass to the limit in the pressure term. 
However let us emphasize at this point that the above a priori bounds do not provide 
any L°° control on the density from below or from above. Indeed, even in dimension 
N = 2, H 1 functions are not necessarily locally bounded. Thus, vacuum patches are 
likely to form in the fluid in spite of the presence of capillary forces, which are expected 
to smooth out the density. It explains why it is so difficult to obtaining the existence 
of global strong solution in dimension N = 2. Indeed in order to prove the existence of 
global weak solution the main difficulty consists in dealing with the quadratic terms in 
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gradient of the density appearing in the capillary tensor. Recently D. Bresch, B. Des- 
jardins and C-K. Lin in [U] got some global weak solutions for the isotherm Korteweg 
model with some specific viscosity coefficients. Indeed, they assume that p,(p) = Cp 
with C > and A(p) = 0. By choosing these specific coefficients they obtain a gain of 
derivatives on the density p where p belongs to L 2 (H 2 ). It is easy at that time with this 
kind of estimate on the density p to get strong convergence on the term of capillarity. 
However a new difficulty appears concerning the loss of information on the gradient of 
u when vacuum existing (indeed the viscosity coefficients are degenerated). It becomes 
involved to pass to the limit in the term p n u n ® u n . That's why the solutions of D. 
Bresch, B. Desjardins and C-K. Lin require some specific test functions which depend on 
the density p. 

In [35], we improve the results of R. Danchin and B. Desjardins in [20] and D. Bresch, 
B. Desjardins an C-K. Lin in [S] by showing the existence of global weak solution with 
small initial data in the energy space for specific choices of the capillary coefficients and 
with general viscosity coefficient. Comparing with the results of [9], we get global weak 
solutions with general test function ip £ C^ D (lR Ar )not depending on the density p. In 
fact we have extracted of the structure of capillarity term a new energy inequality using 
fractionary derivative which allows a gain of derivative on the density p. 
In [25] , we obtain finally the existence of global weak solution for Korteweg system in any 
dimension with large initial data for a specific choice on the viscosity and the capillarity 
coefficients. 

1.1 Mathematical results 

As we explained previously, one of the main difficulty to obtain strong solutions for 
Korteweg's system in very general Besov spaces consists in obtaining regularizing effects 
on the velocity and the density without assuming that the density belongs to L°° . Indeed 
in the general case, we are brought to consider parabolic equations for the velocity and 
for the density with variable coefficients what requires a control on the L°° norms of p 
and on in L°° . To avoid these main restrictions, we will work with a specific choice 
on the capillarity coefficients, namely k(p) = ^ with k > 0. 
In this case, we rewrite the capillarity tensor as follows: 

K(p) =p(VA(lnp) + iv(|VlnH 2 )). 

We refer to the appendix in section [6] for more details on the formal computations on the 
tensor K. We now want to consider the eulerian form of the system (|l.ip when n(p) = — 
with k > 0, we obtain then by dividing by p the momentum equation the following 
system: 

dt In p + u ■ V In p + divu = 0, 
< dtu + u-Vu- ±div(n(p)Du) - ±V(A(p)divu) + VF(p) = VA(lnp) + iVdVlnp] 2 ), 
, (lnp,u) t =o = (In p ,u ). 

(1.4) 

with F(p) is such that — ^ = P ' (p). In the sequel we will use the following definition. 
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Definition 1.2 We set: 

q = In p. 

One can now state the main results of the paper. 

In the first theorem we prove the existence of global strong solution for (jl.4p with small 
initial data when we choose specific viscosity coefficients and pressure terms. We also 
prove the existence of strong solution in finite time with large initial data. The interest 
of this result will be to consider initial data with discontinuous interfaces. 
To do this, we will assume that: 

p(p) = p,p, X(p) = Xp or and with p, X, K > 

It leads to the following system when A = 0: 

dtq + u ■ Vq + divu = 0, 

d t u + u ■ Vu - Au - pVq ■ D(u) - XVqdivu + VF(p) = VAq + iV(|Vg| 2 ), (1.5) 
k (q,u) t =o = (In p ,u ). 



Theorem 1.1 Let N > 2. Assume that p(p) = pp, X(p) = Xp or and P(p) = Kp with 
p,X,K > 0. Furthermore we suppose that: 

AT 1L — 1 

Qo G B p % and u G . 
There exists a time T such that il.5\) has a unique solution (q, u) on (0, T) with: 

Q G (B^) n LUB 2 *Z ), and u G Lf(B 2 ^ ) n 4(^oo )• 

--1 

Furthermore it exists Eq such that if in addition qo G B 2 2 O0 and: 

Ikoll 4-1 4 + ||« || 4-1 < ^o- 

B 2,oo nB 2,oo B 2,oo 

then the solution (q, u) is global and: 

q G L°°(Bf j00 n B£J n L\B*Z n B 2 *Z ), and u G L°°(B 2 % ) n L\B£+ ). 

(1.6) 

Remark 1 We would like to mention that we could easily extend the result of strong 
solution in finite time to the framework of Besov spaces constructed on general LP spaces 
when the initial data are as follows: 

Qo G B p % and u G B^ , 

with 1 < p < 2N . In the case of global strong solution with small initial data in these 
previous spaces, we refer to fT3[ [2S\/ where the ideas could be directly applied. 
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Remark 2 We would like to mention that for the existence of strong solution in finite 
time, we could obtain the same kind of results in non homogeneous Besov space. We now 
want also to emphasize on the fact that when we are working with critical Besov spaces 
for the initial data, we generally are not able to give an estimate of the time of existence 
T. However when we work with subcritical non homogeneous Besov space for the initial 
data, we are able to show that the time of existence for system 81.5\) verifies the following 
inequality: 

T' > °- , 

(1 + IKH + || In poll .f+J' 3 ' 

D P D P 

D p, CO Pi 00 

with 1 < p < 2N and C, (3 depending on the physical coefficients and on the dimension. 
To do this, it is just a matter to show how (ql, u l) <ire small in the critical Besov spaces 
in function of the time and of the initial data. For more details we refer to the proof of 
theorem \1.1[ 

Remark 3 We want to point out that in the theorem \l.l\ we solve the system hl.5\) and 
not the system Ji.^[ j. Indeed in this theorem, we do not assume any control on qo in L°° , 
i.e means that we have no indication on the vacuum. By this way, it is not clear that a 
solution from il.5\) is also a solution from This result proves in a certain way that 

the good variable to consider is not the density p but rather In p. Furthermore we want 
to emphasize on the fact that this result is the first for compressible system where we are 

able to work with such general critical initial data ( in particular uq G -E^oo ) as ^ * s 
case for incompressible Navier Stokes equations (see 111]). In particular it is the first 
time up my knowledge that a compressible system has a velocity which is not Lipschitz. 
Indeed it is generally mandatory to control the L°° norm of the density. 
In some sense, this result means that the Korteweg system with such viscosity and cap- 
illary coefficients is closer to incompressible Navier-Stokes equations than compressible 
Navier Stokes system with constant viscosity coefficients. 

Remark 4 The main interest of this paper is to obtain some solutions which have some 
jump across hypersurface. It means that it could appear some change of phase with 
interfaces of null thickness. Indeed it is remarkable to observe that the density is not 
continuous in your case whereas the interfaces are assumed diffuse. 

Remark 5 We now want to point out the specific choice on the different physical coef- 
ficients. In particular the pressure and the viscosity coefficients have a very particular 
form. Indeed under this form we can check that the variable In p appears everywhere but 
also that there is no non-linear terms in the density such that it would be necessary to 
obtain L°° control on the density for treating these terms. In particular, it is not possible 
to obtain the same type of result for compressible Navier-Stokes equations as we have no 
regularizing effects on the density. Indeed in this case it would be not possible to deal with 

the following term V lnp • in L 1 (5 2 2 00 ). 

In the following theorem we will prove a result of ill-posedness when we consider general 
pressure. Indeed in these case we need to control the L°° norm of the density to deal 
with the pressure term of type P(p) = p 7 with 7 large enough. In the case where 7 is not 
so large, we could extend the theorem \1.1\ by taking advantage of the regularizing effect 
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to estimate Vp 7 in L 1 (B 2oo ). However if we do not assume enough regularity on the 
divergence of the velocity, we loss the L°° control on the norm of the density, it will prove 
in the third theorem. 

We note in the sequel P = (A) _1 Vdiv and Q = (A) _1 curldiv. 

Theorem 1.2 Let N > 2. Let P a suitable smooth function. Assume that p, and A 
verify the conditions of the theorem Furthermore we suppose that: 

JV — — 1 — — 1 

go G B 2 2 x , Fu G B 2 2 2 and Qu G B 2 2 t , 

and that po is bounded away from the vacuum (i.e 3c > such that po > c). There exists 
a time T such that il.5}) has a unique solution (q, u) on (0, T) with: 

q G Lf{B 2 x ) n 4(S 2 y ), P« G Lf(B 2 2 2 ) n L\(B 2 2 ), 

and Qu G LfiB^ ) n L\(B 2 ^ ). 

--1 

Furthermore if in addition go G -B 2 2 i ^en exists Eq such that if: 



\\qo\\ n-! N+\\Fu \\ + Qu 4^ < £o- 
f/ien i/ie solution (g, u) is global and: 

g G L°°(B 2 2 1 n B£) n L 1 ^ 4 " ), Pit G L 00 ^" ) n L\B 2 2 + ), 

and Qu G i 00 ^ ) n L 1 (B 2X ). 

Remark 6 We want to point out here that the solutions of system (1.5)) are also solutions 
of system (1-4^ as q £ L^{L°°) which means that we control the vacuum and the L°° 
norm from the density p. This is then sufficient by composition theorem to prove that 
(p, u) is solution from j 1.1$ . 

We now want to prove that the condition Quo G B 21 is optimal for system (jl.4|) in the 
sense that if we choose a third Besov index which is different from 1 then the problem is 
ill-posed (in the sense of the explosion of norm for the density). 



Theorem 1.3 Let N > 2. Let p and A the viscosity coefficients are as in theorem \l.l\ 
and P(p) = ap. Then there exists a sequel of initial data such that: 

JV N__ -y K — l 

• the sequel (qQ,FuQ,QuQ) is uniformly bounded in (B 2 2 r r\L°°) x B 22 x B 2r with 
1 < r < +oo, 

• we have also: 

jv iY-_i 
q% G B 2 2 v Pu% G Bfa and Q< G B 2 X , 

JV — — 1 — — I 

• (q^ ,Fuq,Quq) tends weakly to (g , Fuq, Quo) in B 2r x B 22 x B 2r , 
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such that the solutions (q n ,u n ) associated to the initial data (go^o) °f theorem \1.2\ verify: 

\\q n (t n , +°° with t n — > 0. 

and (q n ,u n ) converges to a solution (q,u) of U.,5\) . 



Remark 7 This result tell us that the problem is ill-posed for l\1.4\) indeed in this case 
the solution (q, u) of \1.5\i is not a solution of j l.J$ . 



In the specific case where we consider that k = p 2 , we obtain from system (|1 .5|) (when 
A = 0) the following simplified model by assuming that v = u + ^Vlnp (we refer to the 
appendix for more details on the computation or |25j)): 



K 



d t p A/9 = -div(pv), 

» (1.7) 
pdtv + pu ■ Vv — div(pp Vv) + VP(/?) = 0, 

In the sequel, we will set m = pv and q' = p — 1. We have then the following theorem: 

Theorem 1.4 Let k = p 2 . Let P be a suitably smooth function of the density and 

1 < p < +oo. Let mo £ B p p j00 with e > and q Q £ Bp^oo such that po > c > 0. 
There exists then a positive time T such that system \1. 7| ) has a unique solution (q , m) 
with p bounded away from and: 

~ —+e ~ — +2+e ~ — — 1+e ~ — +l+e 

q eC([0,T],B; A )nL 1 T (B p P il ), meC([0,T];B; A ) n L\%T\,B£ X ). 

We assume now that P{p) = ap with a > 0. Furthermore if vq € L°° , Vpo £ L°° and the 
initial data are in the energy space, it means: 

y/p^uo G L 2 , V^p^ £ L 2 and (p - 1) G 

we can i/ien extend the solution beyond (0,T) i/: 

(I _ 1) € Lf{B% +£>l ), (-L - 1) € L^(L l ) and (y/p-l)e ^(L 1 ). (1.8) 



Remark 8 T/ie /irsi condition in il.8\) has to consider as a condition on the high fre- 
quencies and the two last one on the low frequencies. 

In particular, it would be possible to replace the first one, by the following condition: 

(- - 1) € L°°(L 9+£ ) forN = 3 or (- - 1) € L°°(L 3+£ ) for N = 2. 
P P 

Indeed by using energy inequality, Sobolev embedding and by seeing that = — \-VJ~p, 

we can show in this case, that the first one is verified. In some sense these last conditions 
are more practical because they ask only integrability on the vacuum ( and not informations 
in terms of Besov space). However in this case we lost the condition of scaling, we would 
like to explain this point in the following remark. 
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Remark 9 This result is to compare with f32\j where we show a new criterion of blow- 
up for compressible Navier Stokes equations when the viscosity coefficients are constant. 
More precisely we prove that we can extend strong solution if P(p) is in L^(L N+1+£ ). 
In particular we prove that if we have this gain of integrability on the pressure and so 
on the density, we are then able to control the vacuum. In our context here, it is exactly 
the opposite, the criterion of blow up is based on a control of the vacuum i.e — 1) G 
L^(B^ +£ x ). It means that the structure of the viscosity coefficients plays really a crucial 
role in the understanding of the phenomena for compressible fluids. Indeed in our case, we 
can not apply the ideas of f32j because we lost the so-called structure of effective velocity 
introduced in [32]. However as in [32], we are able to transfer the gain of integrability on 
v ( which is due to the fact that we have enough integrability on the pressure ) to a gain of 
regularity on the density. It will be enough to conclude. We can also point out that as in 
]32^ our blow up criterion (where it was the case for the pressure) asks one derivative less 
than the scaling of the equations on (i — 1). It means that in some sense the subscaling 
(one derivative more) required on the initial data (vq G L°° and Vpo G L°° ) is transferred 
in the other sense on the blow-up criterion (indeed we need of one derivative less in terms 
of scaling in our blow-up condition). This fact is very deep. 

Furthermore by following the same arguments than in the proof of theorem \l-4\ it would 
be possible to show a blow-up criterion which would involve ~ in L°° but the initial 
data would be assumed to belong only in Besov space with a very slight subscaling (i.e 

mo G Bp i00 with e > 0, q Q G Bp^ and no more except that the initial data will have 
also to be in energy space). 

Remark 10 We could prove easily that the result of strong solution of theorem ] 1.4\ could 
be extended to critical space in terms of scaling, i.e when we choose e = 0. It will be proved 
in passing in the proof of theorem \l-4\ For a result of global strong solution with small 
initial data for the system \1. 7| ) we refer to \26§ . indeed the proof is a easy application of 
\2b^ excepted that we have to take in account the regularizing effects on the density. 

Remark 11 We could probably extend this previous result for more general pressure 
terms. However it would require more integrability on (- — 1). 

It would be also probably possible to deal with more general viscosity and capillarity but 
there will be certainly an algebraic relation between these quantities. 

Remark 12 In fact it would be easy to improve the previous result concerning the exis- 

. £-1 

tence of strong solution in finite time by working with initial data such that: tuq G B p p 2 

JV 

and q G B p p 2 D L°° such that po > c > 0. To do this we have just to take advantage 
of the first equation in |i. 7\ ) which is a heat equation. We have just to use the kernel 
of the heat equation and write explicitely the solution p to observe that the L°° norm is 
preserved but also the vacuum. 

This article is structured in the following way, first of all we recall in the section [2] some 
definitions an theorems on Besov and Orlicz spaces. Next we will concentrate in the 
section [3] on the proof of theorem 11.11 and 11.21 In section [H we will prove the theorem 11.31 
of ill-posedness. In section [5] we will deal with theorem 11.41 We postpone in appendix 
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(see section [6]) some technical computation on the capillarity tensor and some extensions 
on the previous results. 

2 Orlicz and Besov spaces 
2.1 Orlicz spaces 

We begin by describing the Orlicz space in which we will work: 

L q (R N ) = {/ G Lj oc (R N )/fh\f\<s } e L P (R N ), /l {m > 5} G L q (R N )} 
where S is fixed, S > 0. 

First of all, it is not difficult to check that Lp does not depend on the choice of S > since 
is bounded from above and from below on any interval [Si, 62] with < Si < S2 < +00. 
In particular we deduce that we have: 

f G l\(R n ) if / G L*(R n ) and p, q > e. 

Obviously we get meas{|/| > S} < +00 if / e L^M^) and thus we have the embedding: 

L q p (R N ) C L q p \ (M. N ) if 1 < qi < q < +00, 1 < p < Pl < +00. 

Next, we choose \& a convex function on [0, +00) which is equal (or equivalent) to x p for 
x small and to x q for x large, then we can define the space Lp(R N ) as follows: 

Definition 2.3 We define then the Orlicz space Lp(R N ) as follows: 

L q (R N ) = {/ G L} oc (R N )/*(f) e L 1 (R iV )}. 

We can check that Lp(R ) is a linear vector space. Now we endow Lp(R N ) with a norm 
so that Lp(R N ) is a separable Banach space: 

\\f\\^ m =inf{t>0/ *({)<1}. 

We recall now some useful properties of Orlicz spaces. 
Proposition 2.1 The following properties hold: 

1. Dual space: If p > 1 and q > 1 then (L q {R N ))' = L q ,(R N ) where q = -^ £ , p = 

v 

p-i- 

2. L q = L p + L q ifl<q<p< +00. 

3. Composition: Let F be a continuous function on R such that F(0) = 0, F is 
differentiable at and F(t)\t\~ — > a / at t — >■ +00. T/ien if q> 0, 

F(f)e4(R N ) if feL q (R N ). 

Now we can recall a property on the Orlicz space concerning the inequality of energy. 
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Proposition 2.2 Let P(p) = ap 1 with a > and 7 > 1 f/ien f/ie function H(p) — 11(1) 
is m L 1 (R Ar ) i/ and only i/ (p — 1) € LjJ- 

Proof: On the set {|p — p| < <5}, p is bounded from above, since 7 > 1 we thus deduce 
that jj(p) is equivalent to \p — p\ 2 on the set {\p — p\ < 5}. Next on the set {\p — p\ > 5}, 
we observe that for some v G (0, 1) and C G (1, +00), we have: 

v\p-pV <h{p)<c\ P -pV. 

□ 

2.2 Littlewood-Paley theory and Besov spaces 

Throughout the paper, C stands for a constant whose exact meaning depends on the 
context. The notation A < B means that A < CB. For all Banach space X, we 
denote by C([0,T],X) the set of continuous functions on [0,T] with values in X. For 
p G [1, +00], the notation L p (0, T, X) or Bj,{X) stands for the set of measurable functions 
on (0,T) with values in X such that t ->• ||/(*)IU belongs to L p (0, T). Littlewood-Paley 
decomposition corresponds to a dyadic decomposition of the space in Fourier variables. 
We can use for instance any ip G C°°(R N ) and x £ C°°(M. N ) , supported respectively in 
C = G R^/l < If I < §} and 5(0, |) such that: 

2>(2-'o = i if e/o, 

and: 

Denoting /i = J 7 ^ 1 ^, we then define the dyadic blocks by: 

A t u = ip(2~ l D)u = 2 lN I h(2 l y)u{x - y)dy and S t u = V A k u . 

Formally, one can write that: 

u = ^2 A k u . 

This decomposition is called homogeneous Littlewood-Paley decomposition. Let us ob- 
serve that the above formal equality does not hold in S (R N ) for two reasons: 

1. The right hand-side does not necessarily converge in S' (M N ). 

2. Even if it does, the equality is not always true in S' (M N ) (consider the case of the 
polynomials). 

For the non homogeneous decomposition, we define the dyadic blocks as follows: 
A,u = for I < -2, 
A-iu = x{D)u, 

A,u = ip(2- l D)u = 2 lN [ h(2 l y)u(x - y)dy, for / > 0. 
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Formally, one can write that: 

u = } \ A k u . 



This decomposition is called non homogeneous Littlewood-Paley decomposition. 



2.3 Homogeneous and non homogeneous Besov spaces and first prop- 
erties 

Definition 2.4 We denote by S' h the space of temperate distribution u such that: 

lim SjUj^ +00 = in S . 
Definition 2.5 For s G R, p G [1, +oo], q G [1, +oo], and u G ^'(R-^) we set: 

\\u\\B kq = (^2(2 ls \\A lU \\ LP y)l 

lei 

The homogeneous Besov space B^ q is the set of distribution u in S' h such that < 
+oo. 

Definition 2.6 For s G R, p G [1, +oo], q G [1, +oo], and u G ^'(R^) we set: 

£>'lAH| iP )*)i 



Mr 



P-1 



The non homogeneous Besov space Bp q is the set of temperate distribution u such that 

\\u\\ ,h s < +oo. 

y.q 

In the sequel we will give only properties on the homogeneous Besov spaces but the most 
of then can be generalize for non homogeneous Besov spaces. 

Remark 13 The above definition is a natural generalization of the nonhomogeneous 
Sobolev and Holder spaces: one can show that B^ ^ is the nonhomogeneous Holder 
space C s and that S| 2 * s ^ he nonhomogeneous space H s . 

Proposition 2.3 The following properties holds: 

1. there exists a constant universal C such that: 
C^WuWb^ < ||Vu|| B s-i < C||u||a» ir . 

2. If Pi < P2 and n < r 2 then B s pi Ti ^ J Bp~^ (1/pi_1/p2) . 

' / / 

3. Bp ri B^ r if s > s or if s = s and r\ < r. 



Let now recall a few product laws in Besov spaces coming directly from the paradifferen- 
tial calculus of J-M. Bony (see [5]) and rewrite on a generalized form in [1] by H. Abidi 
and M. Paicu (in this article the results are written in the case of homogeneous sapces 
but it can easily generalize for the nonhomogeneous Besov spaces). 
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Proposition 2.4 We have the following laws of product: 
• For all s € R, (p,r) G [l,+oo] 2 we have: 

\b° i7 . < C(\\u\\ L oo\\v\\ B ^ r + \\v\\L°°\\u\\B° ir ) ■ (2.9) 



uv 



Let (p,p 1 ,p 2 ,r,\i,\ 2 ) e [l,+oo] 2 such that:± < ^ + j^, pi < A 2 , p 2 < A i; | < 

— j- + and - < + We /lave then the following inequalities: 

ifsi + s 2 + iVinf(0, 1 - i - i) > 0, Sl + f < £ and s 2 + § < g then: 

\\uv\\ _ N( j_ + j__ h) < \\u\\ B s lr \\v\\ B s loo , (2.10) 
when si + g = If fresp s 2 + f = we re P lace IMIb^JMIb^ ( res P WAb^J 



N _ N „„j „_ , AT _ TV 

P2 



we tafce 



JMkl,, ^re«P IMIs^nW; */ s i + ~ 2 = n and 82 + AT 
r = 1. 

// Sl + S2 = o, «i € (f - £ , £ - f 1 and -J- + i < 1 iaen: 

J J- ' ^ 7 x v Ai P2 ' PI A2 J Pi P2 — 

\\ UV \\ Nl 1 4- 1 Ii^IMIr"! IMIr"2 • (2.ii) 

If \s\ < — for p > 2 and < s < — else, we have: 
p p p ' 

\\uv\\ B > r < C\\u\\ B ' r \\v\\ k ■ (2.12) 
Remark 14 In the sequel p will be either p\ orp 2 and in this case j = ^ — ^ %fp\ < p 2 , 

Corollary 1 Zei r € [1, +oo], 1 < p < pi < +oo and s such that: 

• »e(-£,£)tf± + £<l, 

• »e(-g+JV(} + £-l),£)tf} + £>l, 

JV 

i/ien we /iave if u £ B s p r and v G B£l t<x n L°° : 

llnw IIrs < C||u||rs lb I 



The study of non stationary PDE's requires space of type L p (0,T, X) for appropriate 
Banach spaces X. In our case, we expect X to be a Besov space, so that it is natural to 
localize the equation through Littlewood-Payley decomposition. But, in doing so, we ob- 
tain bounds in spaces which are not type L p (0, T, X) (except if r = p). We are now going 
to define the spaces of Chemin-Lerner in which we will work, which are a refinement of the 
spaces Lj,{Bp r ). 
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Definition 2.7 Let p G [1, +00], T G [1, +00] and si G R. We sei: 

iiHiz^ (B -) = (E 2 ' rsi ii A ^)ii^(^))"- 

zez 

W^e then define the space L^(B^ r ) as the set of temperate distribution u over (0, T) x R N 
such that \\u\\r Pm s 1 . < +00. 

We set C T (Bp] r ) = Lf(B^ r ) D C([0,T], B^ r ). Let us emphasize that, according to 
Minkowski inequality, we have: 



Remark 15 It is easy to generalize proposition \2.4\ to Lj,{Bp 1 r ) spaces. The indices s\, 
p, r behave just as in the stationary case whereas the time exponent p behaves according 
to Holder inequality. 

In the sequel we will need of composition lemma in L^(Bp r ) spaces. 
Lemma 1 Let s > 0, (p,r) G [1, +00] and u G L p T (B^ r ) n Lf(L°°). 

1. Let F G Wj2 +2 '°°(M jv ) such that F(0) = 0. Then F(u) G L p T (B s p i ). More precisely 
there exists a function C depending only on s, p, r, N and F such that: 

W F ( U )k? T (B kr ) ^ C (\\ U \\L?(L°°))\\u\\ lPT{B s r y 

2. Let F G WQ +3 '°°(R N ) such that F(0) = 0. Then F(u) - F' (0)u G L p T (B^ r ). More 
precisely there exists a function C depending only on s, p, r, N and F such that: 

\\f(u) - F'(o)u\\ ZUBSpr) < c-GMi^oo^inif, {B ,j. 

Now we give some result on the behavior of the Besov spaces via some pseudodifferential 
operator (see [?]). 

Definition 2.8 Let m G R. A smooth function function f : R^ — > R is said to be a S m 
multiplier if for all muti-index a, there exists a constant C a such that: 

v^gr^, \d a f{0\ <c a {\ + \i\) m -\ a \. 



Proposition 2.5 Let m G R and f be a S m multiplier. Then for all s G R and 1 < 

p,r < +00 the operator f(D) is continuous from 5| r to Bp~ m . 

Let us now give some estimates for the heat equation: 
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Proposition 2.6 Let s £ R, (p,r) G [l,+oo] 2 and 1 < p 2 < pi < +00. Assume that 
uq E -Bp r and / € L^ (Bp^ r 2+2 ^ P2 ). Let u be a solution of: 

j d t u - pAu = f 
\ u t=0 = u . 

Then there exists C > depending only on N, p, p\ and pi such that: 

IHIz^+^i) < c{\\u \\ B s r + ^' 1 \\f\\ l P ?{B;7r 2 + 2 /P2) ) • 

If in addition r is finite then u belongs to C([0, T], -B^ r ). 

3 Proof of theorems 11.11 and 11.2 



In this part we are interested in proving the theorems 11.11 and 11.21 of existence of strong 
solutions in critical space for the scaling of the equations. We want to point out that 
the viscosity and the capillarity coefficients are chosen with a very specific form. This 
fact will be crucial in the sequel of the proof in order to obtain estimates on the density 
without assuming a control on the vacuum or on the L°° norm of the density. Indeed 
when fjj(p) = p and n{p) = & with k > then the system has a specific structure and 
provides a new BD entropy (see [25] for more details). 

As a first step, we shall study the linear part of the system (jl.5p about constant reference 
density, that is: 



(iV) 



dtq + divu = F, 

dfU — aAu — bVdivu — cVAq = G, 



3.1 Study of the linearized equation 

We want to prove a priori estimates in Chemin-Lerner spaces for system (N) with the 
following hypotheses on a, b, c, d which are constant: 

< ci < a < Mi < oo, < c 2 < a + b < M 2 < oo and < c 3 < c < M 3 < oo. 

This system has been studied by Danchin and Desjardins in [2D], the following proposition 
uses exactly the same type of arguments. 

— +S ^-1+8 AT 

Proposition 3.7 Let 1 < r < +oo , < s < 1, (go,«o) 6 B^ r x (f? 2 2 r ) N , and 
(F, G) € L\{B*; s ) x {L\{Bi; 1+s )r . 

Let (q,u) e {L\{Bf; s+2 ) n L§?(s| p + ')) x ((^(B 2 f r +s+1 )) JV n (Z§?(S 2 f r + '" 1 ) JV ) &e a 
solution of the system (N), then there exists a universal constant C such that: 
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Proof: 

We are going to show estimates on qi and U[. So we apply to the system the opera- 
tor A; , and we have then: 

d t qi + div-u; = Fi (3.13) 

dm - div(aVni) - V(6divu ; ) - VAg; = G t (3.14) 

Performing integrations by parts and usinf (|3.13)) we have: 



ui ■ V Aqidx = j divui Aqidx, 

d t qi Aqidx + / Fi Aqidx, 
d t Vq r V qi dx+ [ VFfVqidx, 

JR N 

I \Vqi\ 2 dx- / Vqi.VFdx. 

JR N JR N 



ld_ 
2~dt 



Next, we take the inner product of (|3.14|) with U[ and by using the previous equality, we 
have then: 

lll(\\ u l\\h+ [ \Vqi\ 2 dx)+ [ (a|Vu,| 2 + 6|divn i | 2 )(ix = 
2 dt Ji" Jrn 

G l .u l dx+ [ Vqi.VFidx. ( 3A5 ^ 



In order to recover some terms in Aqi we take the inner product of the gradient of (|3,13D 
with ui, the inner product scalar of (|3.14|) with and we sum, we obtain then: 

^ / Vq l .u l dx+ / {Aqifdx=f (Gj.V« + |divn,| 2 + Ul .VF t 

dt Ju N Jr n Jr n (3.16) 

— aVui : V 2 qi — bAqid\vui)dx. 

Let a > small enough. We define: 

kf = \\uif L 2 + / (K\Vqi\ 2 + 2aVqi.ui)dx . (3.17) 
Jr n 

By using (|3.15p . ()3.16p and the Young inequalities, we have by summing and the fact 
that a is small enough: 

1 d 1 f 

oT+rf + o / (a\V Ul \ 2 +ab\Aq l \ 2 )dx<\\G l \\ L 2(a\\Vq l \\ L 2 + \\u l \\ L 2) 

2dt 2 J r n (3.18) 

+ l|V.F}|| L 2(a||u/|| L 2 + ||V^|| L 2). 

For small enough a, we have according (|3.17p : 

\k 2 < \\ui\\ 2 + [ R\V gi \ 2 dx < hf . (3.19) 
2 J r jv 2 
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Hence according to (|3.18p and (j3. 191) : 

\jk} + K2 2l k? <Ch + ||VF ; || L2 ). 

By integrating with respect to the time, we obtain: 

hit) <e- K22l %{Q) + C f e-^ 2i ^)(||VF,(T)|| L2 + {{G^H^dr . 

Jo 

After convolution inequalities imply that: 

\\k\\ L p([o,T]) <C{2- 2 ik l (0) + 2- 2l{1+l P -n ) \\(VF l ,G l )\\ LPTl(L2) ). (3.20) 
Moreover we have: 

C~ x ki < \\Vqi\\ L 2 + \\ui\\ L 2 < Ck h 
Finally multiplying by 2 V 2 p , taking the V norm and using ()3.19j) . we end up with: 

||(Vg,n)|| n 1+s+ 2 < ||(VF,G)L ^_ 3+s+ ^ + ||(V<fo, «o)|| 

It conclude the proof of the proposition. □ 
We now want extend the result of the proposition 13.71 to the case where we include the 
pressure term inside of the linearized system. This point will be crucial in order to deal 
with the existence of global strong solution with small initial data. Indeed in this case it 
is very important to take in account the behavior in low frequencies and so the pressure 
term which is local. More precely we will consider the following linear system: 



(iVl) 



dtq + divu = F, 

dtu — aAu — 6Vdivu — cVAq + dVq = G, 



We now want to prove a priori estimates in Chemin-Lerner spaces for system (-/VI) with 
the following hypotheses on a, b, c, d which are constant: 

< ci < a < Mi < co, < c 2 < a + b < M 2 < co, < c 3 < c < M 3 < co 

and < C4 < d < M4 < 00. 

This system has been studied by Danchin and Desjardins in [20J, the following proposition 
uses exactly the same type of arguments. 



1+s 

2,r 



Proposition 3.8 Let 1 < r < +00 , < s < 1, (q ,u Q ) G (B£ r D B, 2 
(B*- 1+S ) N , and (F, G) G L* (B£ 1+8 n B^*) x (L^B^ 1+S )) N . 

Let ( q ,u) e {LUB 2 2 ir + 8+1 n B*; s+2 ) n L??(^ r r 1+S n i^ r + )) x ((L^ r +S+1 ))^ n 



AT 



(L^ (B 2 2 r +S 1 ) Ar ) ^ e a solution of the system (-/VI), i/ien i/iere exists a universal constant 
C such that: 

ll(v ^'" )ll r.(4-- ) nr P(B J-, £ c (IK^o, 9 c,«o)ll 4+ . + IW'^l 

Proof: 

It suffices to follow exactly the same lines as the proof of proposition 13.71 
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3.2 Proof of the theorem ll.lt 

We now are going to prove the existence of strong solutions in critical space for system 
(jl.5p . In particular we recall that the main interest of theorem 11.11 is to allow discontin- 
uous initial data for the density, such that we can authorize discontinuous interfaces. 

Existence of solutions 

We use a standard scheme: 

1. We will use a classical iterative scheme to constructed a sequence of approximated 
solutions (q n ,u n ) on a bounded interval [0,T] which depend not on n. We will get 
uniform estimates on (q n ,u n ) in: 

£ T = (c^Jn4(< )) x (C T (B 2 % ) n )). 

2. We will prove that the sequence (q n ,u n ) is of Cauchy and converges to a solution 

of flr 



First step 

We smooth out the data as follows: 

<?0 = S nQo, Uq = S n u and f n = S n f. 

Note that we have: 

V/ G Z, ||A^o < ||A;<7o||l2 and \\q%\\ n < \\q \\ n , 

and similar properties for Uq and f n , a fact which will be used repeatedly during the 
next steps. Now, according [33], one can solve (|1.5p with the smooth data (qQ,UQ,f n ). 
We get a solution (q n ,u n ) on a non trivial time interval [0,T n ] such that: 

fe^T^B^ni^B^ ) and u n e C([0,T n ), )n4„(B^ + ). (3-21) 

Uniform bounds 

Let: 

q n = QL + q n , u n = u L + u n , 
where {qiiUi) stands for the solution of: 

\d t q L + divu L = 0, 
\d t u L - Au L - nV{Aq L ) = 0, 

supplemented with initial data: 

<?l(0) = go , ul(0) = u . 



(3.22) 
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Using the proposition 13.71 we obtain the following estimates on (ql,ul) for ah T > 0: 

gL eC([0,T],B 2 ^)n4(< ) and u L 6 C([0, T], )n4(B£ ). 

We let (g°,n°) = (0,0). We now want study the behavior of (q n ,u n ) where (q n ,u n ) are 
the solution of the following system: 



(Ni) 



' d t <f + div{u n ) = F n _i, 
< <9 t u n - Au n - KV(A<f ) = G n _i, 
k (q n ,u n )t=o = (0,0), 

where: 

F n _x = - u^ 1 • Vq n -\ 

= -u L -Vq L - u n ~ l ■ Vq L - u L • Vg™" 1 - u™" 1 • Vg™" 1 , 

G n _! = - {u^y.Vu 12 - 1 + /iVg™- 1 • Du n ~ l + AVg™" 1 divu™- 1 + IvdV^" 1 ) 2 ) - KVq n ~ l . 
1) First Step , Uniform Bound 

Let e be a small parameter and choose T small enough so that by using the estimate of 
proposition 13.71 we have: 



We are going to show by induction that: 

(V n ), \\(q n ,u n )\\ FT <V~e. 
for e small enough with: 

F T =(C([o,r],B 2 yn4(< )) x {C([0,T],B 2 % )n4(< ))• 

As (g ,^ ) = (0,0) the result is true for n = 0. We now suppose {V n ~i) (with n > 1) 
true and we are going to show (V n )- Applying proposition 13.71 we have: 

\\(q n ,u n )\\ FT <C\\(VF n ^,G n ^)\\ Ti n w (3.23) 

Bounding the right-hand side may be done by applying proposition 12.41 lemma [T] and 
corollary [TJ We begin with treating the case of ||i^ n _i|| ri ,„jv/2 % , we have then: 

\\UL • ^ Ql\\ fl m N/2s < \\ul\\j 1 rR JV/2+lJ gi / |Ur [x> , R iV/2, + | 9i[ 71 , R N/2+2,J \Ul\ 7^ , R JV/2-l v 



Similarly we obtain: 



UL-Vq n 1 \\fl m N/2, < \\ulWj! m N/2+lA\q n 1 ||roo/ R JV/2,+||g n || 4 JV/2+3 JkilU. , iV/2-Il 
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Il^-Vgzll^! (R N/2. < \\u n 1 \\ 4 JV/2+1 JQl\\ m/2+1 +ll?i||fl (R N/2+2 ) \\u' n 1 1 1 ~ N/2-1 

and: 

ll^-rn — 1 \7~n— 111 I Ill ||=ra— In i II ^n— In 1 1 ^rj, — In 

F " V 9 Tl/ R iV/2^S « \\f llR N/2+l. \\q 7oo rR JV/2^+ 9 7l /R^/2+2j « Too/r^-In- 

L T\ B 2, oa ) T v 2,oo J L t( B 2,co) L t( B 2,oo ) T \ 2,oo J 

By using the previous inequalities and ('He), we obtain that: 

11^11^(^/2) < C(2^ e + 2ei + 2^ie\ + 2e). (3.24) 
Next we want to control ||G n ||_ jv_ 1 . According to propositions 12.41 corollary [1] and 



13. 1\ we have: 



f„, n— 1^* V7», n — 111 <* IL, ra — 111 IL, n— In 

It .Vtl H , >, 11 4 1,1 U N 1 



T^ B 2,oo ) L t( B 2^oo^) L t( B 2^oo^) 

n-l|2Mi <- II lv7„n-l|2|i 



HVdVg^niL <|||V<f 



A" 



L t( B 2,oo > L T\ B 2,vo 



111 ~ 4 N -I- 3 111 \\~ ^-t- 1 • 

We proceed similarly for the other terms and we obtain by using (|3.23|) and the different 
previous inequalities: 

\\(q n+1 ,u n+1 )\\ FT < Cy/e(y/eAo + ^/e + e^). 

By taking T and e small enough we have (V n +i), so we have shown by induction that 
(q n ,u n ) is bounded in Ft- 

In the case where we want to obtain global strong solution with small initial data, we 
need to take in account the low frequencies. That is why we will include the pressure 
term in the linear part and we will use proposition 13.81 to obtain estimate on (qx, ul) 
where (<7l,ul) here stands for the solution of: 

(d t q L + divu L = 0, 

[dtUL - Au L + KVq L - KV(Aq L ) = 0, 

In particular it explains why we need of additional regularity on the initial density, i.e 

qo G ^200 • The res t °f the proof for this case is similar to the case with large initial 
data. 

Second Step: Convergence of the sequence 

We will show that {q n , u n ) is a Cauchy sequence in the Banach space Ft, hence converges 

to some (q, u) £ Ft- 

Let: 

5q n = q n+l - q n , 5u n = u n+1 - u n . 
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The system verified by (8q n ,5u n ) reads: 

f d t 5q n + div5u n = F n - F n _i, 

d t 5u n - fiA5u n - (A + fi)Vdiv5u n - VA5q n = G n - G n _i, 
^ tf g »(0) = , <5u n (0) = 0, 

Applying propositions 13.71 and using (V n ), we get: 

\\(6q n , Su n , 5T n )\\F T < C(\\F n - F n _a|| N/2 + ||G n - G^W (B w-u). 

And by the same type of estimates as before, we get: 

U5q n ,6u n )\\F T <C^-e{l + A f\\{8q n -\8u n - l )\\ FT . 

So in taking e enough small we have that (q n ,u n ) is Cauchy sequence, so the limit (q,u) 
is in Ft and we verify easily that this is a solution of the system. 

Third step: Uniqueness 

Now, we are going to prove the uniqueness of the solution in the following space: 

F- = (C([0,T],B 2 %)n4(< )) x {C([0,T],B 2 % )nL 2 T (B 2 %)). 

Suppose that (qi,u%) and (q2,u 2 ) are solutions with the same initial conditions, and 
(<7i,ui) corresponds to the previous solution. 
We set then: 

6q = q 2 — qi and Su = u 2 — U\. 
(5q, 5u) satisfy the following system: 

' dt5q + div<5u = F 2 — F\, 

< d t Su - fiA5u - (A + /j,)Vdivdu - VA5q = G 1 -G 2 , 
^ 5q(0) = , 5u(0) = 0. 

We now apply proposition 13.71 to the previous system, and by using the same type of 
estimates than in the part on the contraction, we obtain: 



\\(Sq,6u)\\M < (\\qi\\ n +1 n +1 , tf+i +||«2||_ 9/ n ) 

x \\(5q,6u,ST)\\M. 

We have then for T% small enough: (&/, 6u) = (0, 0) on [0, T\\ and by connectivity we 
finally conclude that: 

Qi = <?2, u\ = u 2 on [0,T]. 

□ 
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3.3 Proof of the theorem 11.21 



In this section we have to deal with general pressure. It means that we need of addi- 
tional information to control the pressure in Besov space. Indeed in order to use some 
composition law of paraproduct on the pressure it is necessary to control the L°° norm 
of the density p. To do this we need of additional regularity assumption on the initial 
data to be able to estimate the L°° norm on the density p. 

The proof follows the same lines as the proof of theorem 11.11 except that we obtain new 
estimates on (ql, ul), indeed by using proposition 13.71 (ql^l) is in Et with: 

~ M. ~ iv ,p ~ iL—-\ ~ Ji-n 

Sr = (c([o,n^)n4(< )) x [C([0,T\,B& ) n L\{B£ )). 

However in order to estimate via the mass equation the density in L°° , we need to 

control dfvu in L 1 (L°°).That is why as Quo belongs to B 2 2 1 , we are able to prove that 
~ u. 

divuL £ L 1 (B 2 2 1 ). Indeed applying the operator div to the momentum equation of system 
(|3.29p and letting vl = divui, we obtain the following system: 

d t c + Av L = 0, 

d t v L - (2p + X)Av L - k,Ac = 0, (3.26) 
(c,v L ) /t=0 = (Alnp ,divn ), 

with c = A In p. We now recall a lemma obtained in [20J . 

Lemma 1 Let sGl and r G [l,+oo]. Suppose that vq G -B|r- Then the system (??) 
has a unique solution vl in C(i?| r ) Pi L 1 (i?2 and we have: 

KHz~(bs iP ) + II^IIli(b=+ 2 ) ^ c W v o\\Bi r , 
where C depends on p, A and k. 

As a consequence of lemma HJ we obtain that divu^ belongs to L 1 (i3 2 2 i). Unfortunately 
it is not sufficient to hope controlling the density in norm L°° . Indeed via the transport 
equation, we need also to control the divergence of the velocity in L 1 (L°°). To do this, it 

suffices to estimate u n in L 1 (S 2 2 1 ). The idea is then to prove that u n is more regular on 
the third index than ul- This type of result is well-known in the case of Navier-Stokes 
(see [12j). we now want to adapt the spirit of this result to our case. 
We now have to prove that (q n ,u ?l ) is bounded in E)p with: 

4=(C([0,T],B 2 yn4(< )) x {c([o,t],b 2 >- )n4(< ))■ 

It means in particular that we have a type of regularizing effects on {q rl ,u n ) on the 

third index of the Besov spaces. More precisely we are going to show by induction as in 
theorem 11.11 that: 

(Vn), \\(q n ,u n )\\ E i T <e. 
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for e small enough. 

As (if, u°) = (0, 0) the result is true for n = O.We now suppose (V n -i) (with n > 1) true 
and we are going to show (V n )- Applying proposition 13.71 we have: 

\\(q n ,u n )\\F T <C\\(VF n ,G n )\\ Ti n (3.27) 

Bounding the right-hand side may be done by applying proposition 12.41 lemma [TJ and 
corollary [TJ We want only to give an example of how to control ||G n |L .v , • We will 

^Woo ) 

treat the case of ■u n ~ 1 )*.Vu™ -1 and V(\Vq n ^ 1 \ 2 ). According to propositions 12.41 corollary 
[TJ and 13.71 we have: 

,n—l\* v7„ vt, — 1 1 1 ll„.fi— 1|| ll„ ft— 1 

■ VU N , > 



lit j .VU N >- \\U N _, \\u 

M ) L T ( B 2% ) "^ftYl 



livdv^- 1 ! 2 )!! ^ < lliv^- 1 ! 2 !! * , 

L T\ B 2,oa ) L T\ B 2,l) 

< \\Vq n - l \\ 4 iv.i HVg"" 1 !! iv i , 
<" 1 1 — 1 1 1 II „rt- — • 111 

||9 ||~4 iV + 3 ||9 ||_ N + l . 

Lt( B 2,2 ) L T {B 22 ) 

By following the same lines than the proof of theorem 11.11 we can easily conclude that 
(V n ) is verified. The only difficulty consists in treating the non linear term coming from 
the pressure. To deal with this term, we need to prove that p n belongs uniformly in L°°. 
From lemma [T]and proposition 13.71 we obtain that: 

ql G L°°(B£j) n L (B^i ), div< G 2 ) n L 1 ^) 

and u£ G L°°[B% X ) n ^(^l )■ 

It means in particular as u n is in L l (B^ ), and from f[3T28l) we conclude that divu n is in 
~ — +i 

L 1 (i? 2 2 i )• It is then sufficient to prove that q n is in Lj? {L°°). Indeed via the transport 
equation, we have: 

1 1 9*1 1 1,°° < lko|U°° + ||divit n || L i( L c=o). 
We can the deal with the pressure term as follows by using proposition [TJ 

„ 1 



V 

T N 



P L t( B 2,1 ) L t( B 2,\) 

We can now conclude easily the proof of theorem 11.21 by following the same lines as in 
theorem 11.11 

In the case of the existence of global strong solution with small initial data, we need to 
take in account the low frequencies. That is why we will include the pressure term in the 
linear part and we will use proposition 13.81 to obtain estimate on (9l,u_l) where (qL,uz) 
here stands for the solution of: 

(d t q L + divu L = 0, 

1 d t u L - Au L + KVq L - KV(Aq L ) = 0, 
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In particular it explains why we need of additional regularity on the initial density, i.e 
--1 

qo G ■ The rest of the proof for this case is similar to the case with large initial 

data. □ 



4 Proof of theorem 11.3 



We now want to prove theorem 11.31 which is an ill-posedness theorem, in the sense that 
we have an explosion of the L°° norm of the density in a arbitrary small time. To do this 
we have just to choose a sequence of initial data (c/o,Uo)neN such that: 

K — i £L — — 1 

• the sequel (q^ , Pkq , Q^o) is uniformly bounded in (i? 2 2 r n -B 2 2 r fl L°°) x £? 2 2 2 x 

B 2 \ with 1 < r < 2, 



we have also: 



7 £ G B* v P< G B£ and Qu% G Bfa , 



• ($,Pt#,Qit&) tends weakly to (q ,¥u ,Qu ) which is in (B 2 2 r n B£ r n L°°) x 

D 2 x R 2 

"2,2 x D 2,r ' 

Furthermore we ask that there exists a sequel (tn)neN such that t n — > n ^+ 00 and that 
the solution (q^,^) of system (|3.29p with initial data (<7o,iig) verify: 

||divu£||ij fL oo) -^n^+oo +oo. 

By theorem 11.11 we can show easily that it exists T > such that it exists strong solution 
(q n , u n ) verifying (|1.6j) with initial data (qq , Ug) on (0, T). Furthermore by using the same 
arguments than theorem II. 2\ we can show that: 

q n = ql + q a and u n = u n L + u n , 

such that (q n ,u n ) is uniformly bounded in the following spaces: 

G L?^) n L^^ 2 ), G L°°(B£ X ) n L (B^x )■ (4.30) 

and such that by using proposition I3.8I {(fi^vlf) is uniformly bounded in the following 
spaces: 

g£ G L™m r ) n L\B* +2 ),ul G L5?(SA _1 ) fl )• (4-31) 



We now want to consider the transport equation of system (jl.5p . we have then: 

dtq n = -divu™ - u n ■ Vq n . 

We deduce that: 

q n (t,x) = - / divu n (s, x)ds - / u n ■ Vq n (s, x)ds. 
Jo Jo 
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By paraproduct we can show that: 




N 



We deduce that u n ■ Vg n is uniformly bounded in L 1 (B 2 2 1 ) and then : 




Jo 



N 



As divu n = divu n + divu^ with divu™ uniformly bounded in L^B^i) Ly(L°°). 



Furthermore we recall that: 



l|div<|| L : 



n— >+oo ~l~00. 



It finally shows that: 



n- 



OG 



+ CO 



which is the desired result. 



□ 



5 Proof of theorem 11.41 

For the proof of the existence of strong solution in finite time on (0, T), it is a basic 
application of theorem 5 in [20] except that we do not ask condition of smallness on the 
initial data. To do that we proceed as in [33]. We just want to mention that we consider 
the momentum variable in order to choose some p arbitrary big. Indeed it is crucial when 
we want to deal with the term div(um) that the quantity um be inside of the divergence. 
It allows us to apply paraproduct laws. 

We now want to prove the more interesting part of theorem 11.44 it means the blow-up 
criterion. The main idea is to obtain a gain of integrability on our solution (p,v). By 
passing we can observe that (In p, u) is also a solution of system ([1.50 because we have a 
control on In p. After that we will use this gain of integrability on v in order to obtain 
additional regularity on the density. It will be enough to conclude. 

In the sequel we just will consider the case N = 3, the case N = 2 follows the same lines. 
Gain of integrability on v 

We multiply the the momentum equation of (|1.7p by u|-u| p-2 and integrate over R , we 
obtain then: 




(5.32) 



Next we observe that: 




26 



We get then as div(pu) = —dtp and by using (|5.32|) : 

- / d t (p\v\ p )dx + [ p\v\ p - 2 \Vv\ 2 dx + {p-2) [ p(V di(\v\ 2 ) 2 \v\ p - A dx 

P Jr n Jr n Jr n i ^ ^ 

+ f \v\ p - 2 v -Vp^dx = 0. 
Jr n 

We have then by integrating over (0, t) with < t < T: 

- f (p\v\ p ){t,x)dx+ f f p\v\ p - 2 \Vv\ 2 (t,x)dtdx 
P Jrn Jq J r n 



+ 



(p-2) f [ p{Y j d l {\v\ 2 ) 2 \v\ p -\t,x)dtdx< - I (p \vo\ p )(x)dx (5.34) 
Jo Jrn s-f p J r n 



+ 1 I I \v\ p ~ 2 v-Vp J (t,x)dtdx\. 
Jo Jr n 



By integration by part we have: 
ft 

la Jr n 



v\ p ~ 2 v ■ Vp 7 (i, x)dtdx = - [ [ dw(\v\ p - 2 v)ap(t, x)dtdx. 

Jo Jr n 



We have then: 

div{\v\ p - 2 v) = \v\ p - 2 div(v) + (p - 2)\v\ p - A v ■ (v ■ Vv). 

In the sequel we will note p = N + e with e > and we will use the fact that as ( 1) 

is in Lf(B° N+£A (R N )), it is also in Lf(L N+£ (R N ). We can prove that (i - 1) is in 
L°°(L 2p ) because we have: 

(^=-l) 2p =(( i -l) + 2(l-^=)) P , 
y/P P y/P 

and: 

||-^-l|U*<ll(--l) + 2(l-4p||i, 

< c(i +\\{- p - i)iii +eii(i - ~^)\\ L 2 P +\\-j=- in!)- 



We now recall that by (|5.34p . we control ^/pV(|ti| 2 ) in L^(L 2 ) and as — 1) belongs 

I = I + JL 

q 2 T 2p- 



to L°°(L 2p ) we can show that V(|u|§) is in L^(L 9 ) with \ = \ + 1 



By Sobolev embedding we show easily that: 

HMi- 1 !! * <c||v(|^|l)|U L2) , 

and similarly that: 



I^ILj^^^ < c||v(Ht)|U iq) . 
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with i = I + ^ — By energy inequality ([1.30 Vy/p belongs to L^(L 2 ). By Sobolev 
embedding, we deduce as (yfp — 1) is in L°°(L 1 ) that (y^o- 1) is in L°°(L 6 ). 

Finally, we have proves that y / p(y / p|f| 2 _1 Vu)|w| 2 _1 is in L 1 (L a (M Ar ) with a > 1. We 
recall that: 

V^VpM^V^M* -1 = (VP" 1)(\/p|w|2 _1 Vu)|u|2 _1 + (Vp|w|2 _1 Vu)|u|2 _1 . 

To do this we just apply Holder's inequalities ,as by ([5.340 we have a control on yfp\v\ 2 _1 Vi> 

mLl(L 2 (R N )), on (y/p-1) mLf(L 6 (R N ) and |-^| f ~ 1 mL^(L^). We have then the 
following condition: 

1 p-2 1 1 p-2,1 1 1. 1 1 

- = - + - + - = - + ) + ~ + - < 1, 

a pqi 2 6 p y 2 2p 3 J 2 6 ~ ' 

because p = 3+e > 3. By Sobolev embedding, we can exactly prove that y/p(y/p\v \ 2 _1 Vu)|w| 2 _1 
is in fact in L^L 1 ) with: 

/ \v\ p ~ 2 v ■ VP(p)dxdt\ < \W\v\S)\fef m {l + ||- - 1|U» (LP) 
o Jk^ v ; P 

+ II — 7= - illi 00 ^ 1 ) + IIvp - 1 ll J L°°(i: i )) /3 7 

VP 

i 

with /3 big enough. By bootstrap in ()5.34p we obtain finally that ppv is in L°°(L P ) for 
any 1 < p < +oo. 

Blow-up criterion 

i 

We have then obtain that than for any 1 < p < +oo, p?v belongs to L°°(L P ). We now 
want to consider the first equation of ([1.70 : 

dtq — —Aq = —div(pv), (5.35) 
P 

with q = p — 1. We would like to take advantage of the gain of integrability that we have 
obtained on the velocity v in order to show regularizing effect on the density. The goal 
is to transfer the information on the integrability of v (which is a subscaling estimate) 
on the density p. More precisely we have by proposition 12.61 for any 1 < p < +oo: 

lls'lll-^) < Cdkollloo^) + IHIl^BO^)). (5.36) 

We know by energy inequalities ([1.30 that: 

Go- 1) G L°°\Ll) and V^p € L°°(L 2 ). (5.37) 

We then split the product pv as follows: 

i i ii /ii 

pv = (p p — l)ppy + ppv = f(q )ppy + ppy, 
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with f(x) = (1 + x) p — 1. We now need to prove that pv is in L^?(i?p j00 ). As p?v is 
in L^?(L P ), we have then: 



(5.38) 



11/(9 )P p v\\l™(Lp) < \\p p v\\L™(LP)\\f(q 

<C||P p v||L5?(LP)(l + ll? lli/(L-))> 

But by interpolation as p is big enough, we have for : 

Ik'lU^-) < lis ll^^jll? lISp^). ( 5 - 39 ) 

with < a < 1. By injecting ()5.38|) and (|5.39|) in ()5.36|) . we obtain: 

Us ll^?^.) ^ ^(llflbllr-^,.) + Wp'Awwi 1 + H« ll^ iao )H« II2?(£»))- ( 5 - 4 °) 

As by (|1.3p and Sobolev embedding, H^'H^oo^aj is finite, we deduce from (|5,40p and from 
Young's inequalities that q is in L^i^Bp^) for any N < p < +oo. It means that we have 
obtain that Vp G L^?(Sp i00 ). As in theorem 1 1.4} we assume that p G L^'(B < ^ +£ x ) with 
e > 0, we have by paraproduct law (see proposition I2.4p when we are the critical case 
where the sum of the regularity index in nul: 

||Vmp|| Zoc( ^ 2i0c) < C(l + ||- - n~ L ^BO N+£tl) )\\^Ph^B0^ 

with ~ = | + TvTe = tv+- (when we work with p large enough). We obtain then that 
Vlnp G L°°(B° N+ e J and that lnp G ^(B^ J. 



Furthermore we have obtained that ppy belongs to L°°(L P ) for any 1 < p < +oo and 
that Vlnp G Lf{B° N+ ^ oo ). As u = v - ^Vlnp and the fact that (± - 1) G L^(B% +£ ), 

we obtain easily by following the same lines that for p large enough u is in L^(B^ + e OQ - 
To summarize what we have obtained, we have: 

ueL^(B° N+hoo ) and lnp G L§?(*jv+S (0O ). (5.41) 

We recall that the solution (p, u) of system (|1.7p on (0, T) are also solutions of system 
p.4p in the sense that (q,u) is solution on (0,T) of system (jl.4p . 

Easily we can prove that if no G B^ + ^ oo ) and lnpo = Qo G ^n+s. oo then the system 
(|1.4|) has a strong solution (g, u) on (0, T ) with: 

C 

T > 



u + ikiibo + ii lnpoiis 1 F y 



This is an easy consequence of the fact that the initial data are choose subcritical, indeed 
in this case we have just to use proposition [321 to prove that ml) are small in function 
of the time T in the critical norm. 

It means that there exists a time T > c > 0, where c depends only on the physical 
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coefficients and of subcritical initial data. We can construct by theorem 11.11 a solution 
(qi,u±) on (T — a,T — a + T ) with initial data (q(T — a),u(T — a)) (here a < T ) which 
verifies (|l,6p (in fact with a bit more regularity as we are subcritical) with Lebesgue 
index p = N + |. The only difficulty is to prove that on (T — a, T) we have: 

0?i, ui) = i,q,u). 

To do this, it suffices only to use the uniqueness part of theorem 11.11 It concludes the 
proof of theorem 11.31 □ 

6 Appendix 

In this appendix, we just give a technical lemma on the computation of the capillarity 
tensor. 

Lemma 2 When n(p) = — with n > then: 

divK = K P (VAlnp+ ^VGVlnpl 2 )). 

and: 

divK = Kdiv(/)VV In p). 

Proof: We recall that: 

divK = V(pK(p)Ap + -{k( P ) + pk (p))\V p\ 2 ) - div(«;(p)V/) ® Vp). 
When k(p) = ^, we have: 

divK = kVA/9- Kdiv(ivptg) Vp). (6.42) 

P 



But as: ^ 

Ap = pAlnp+ -|Vp| 2 , 
P 



we have by injecting this expression in (|6.42p : 

divK = npVAlnp + KVpAlnp + KV(-\Vp\ 2 ) - Kdiv(-Vp®Vp). (6.43) 

p p 

As we have: 

1 Ik 

fixhv(-Vp® Vp) = KAlnpVp + V(-|Vlnp| 2 ) pV(|Vlnp| 2 ). 

p p 2 

It concludes the first part of the lemma. 
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We now want to prove that we can rewrite (??) under the form of a viscosity tensor. To 
see this, we have: 

div(pV(Vlnp)) 3 - = ^2 di(pdij In p), 

i 

= ^2 [dipdij In p + pduj In p] , 

i 

= p( A V In p)j + ^2 P®i ln P 9 3 9 i ln P) ' 

i 

= / >(AVln^ + ^(V(|Vln / >| 2 ))„ 
= divK. 

We have then: 

divK = Kdiv(/9VVlnp) = «div(p£)(V lnp)). 

□ 
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